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We present a comprehensive analysis of the halo model of cosmological large to small-scale struc- 
ture statistics in the case of warm dark matter (WDM) structure formation scenarios. We include 
the effects of WDM on the linear matter power spectrum, halo density profile, halo concentration re- 
lation, halo mass function, subhalo density profile, subhalo mass function and biasing of the smooth 
dark matter component. As expected, we find large differences at the smallest physical scales in 
the nonlinear matter power spectrum predicted in the halo model between WDM and cold dark 
matter even for reasonably high-scale WDM particle masses. We find that significant effects are 
contributed from the alteration of the halo density profile and concentration, as well as the halo 
mass function. We further find that the effects of WDM on the subhalo population are important 
but sub-dominant. Clustering effects of the biasing of the smooth component in WDM is not largely 
significant. 
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I. INTRODUCTION 

The theory of large scale structure growth and for- 
mation from Gaussian, adiabatic initial conditions is in 
remarkable agreement with observed large scale struc- 
ture in the universe at scales above the Galactic scale. 
With the success of the Wilkinson Microwave Anisotropy 
Probe p] at determining the amplitude and shape of the 
matter power spectrum, and upcoming results from the 
Planck Probe, the primordial linear clustering of mat- 
ter will be known to better than 1.1% in the amplitude 
and 0.22% in the slope of the power spectrum. When 
combined with galaxy surveys such as the Sloan Digi- 
tal Sky Survey (SDSS) and 2-degree Field, the precision 
can be even further enhanced [2] . These precise determi- 
nations of the clustering power spectrum, coupled with 
the ansatz of cold dark matter (CDM) with a primordial 
power law spectrum leads to a precise prediction for the 
power spectrum to arbitrary small scales of relevance to 
galaxy formation. 

However, the small-scale clustering of dark matter can 
be largely different at the small scales inaccessible to 
CMB measurements as well as galaxy surveys, yet are sig- 
nificant in the role of galaxy formation. The small-scale 
clustering of matter can be suppressed by thermal mo- 
tion of the dark matter, i.e., warm dark matter (WDM), 
allowing the dark matter to free stream out of primordial 
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potential wells. It can also be altered in extended infla- 
tionary models [3] , and in the case of charged or neutral 
decaying [4] dark matter progenitor particles. 

The suppression of power on small-scales has been pro- 
posed as a solution to a number of problems in galaxy for- 
mation: first, the reduction of satellite galaxy halos [5], 
second, the reduction of galaxies in voids [6] [7], third, the 
low concentrations of dark matter in galaxies [8j , fourth, 
the angular momentum problem of galaxy formation [5] , 
fifth, the formation of disk-dominated galaxies [TO] . 

The discovery of a number of dwarf galaxies with the 
SDSS, and their confirmation as harboring massive dark 
matter halos [UJ, alleviates the first problem, and even 
more, constrains the excessive suppression present in 
WDM models 12J. In addition, the small-scale power 
spectrum of dark matter inferred from the Lyman-a for- 
est toward distant quasars is consistent with CDM down 
to sub-dwarf galaxy scales [13] , though the thermal his- 
tory of the intergalactic medium consistent with these 
findings is difficult to reconcile with other measures [14] . 
Overall, however, there is considerable interest in WDM 
as a potential structure formation scenario, and as a 
method to test the CDM ansatz. 

In this paper, we construct a complete analytic model 
for the clustering of dark matter in two-point statistical 
measures based on the halo model of large scale struc- 
ture. For a review of the halo model, see Ref. [U]. The 
forms developed here may be extended to higher order 
statistics as well. The halo model has been shown to 
be an effective and accurate model of dark matter clus- 
tering built upon analytic and numerical structure for- 
mation statistics. And, it can be extended to models of 
galaxy [IB], gas [T7] and galaxy velocity statistics [T8] . 
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It has been shown to be consistent with broad classes of 
galaxy clustering measures over a wide range of redshifts, 
c.f. Ref. QJ]. 

The halo model has been employed to investigate the 
effect of baryons [50] and massive active neutrinos [5T] 
on weak lensing statistics. Recently, the effects of WDM 
in weak lensing has been studied through estimates of 
WDM's alterations of limited components of the halo 
model. In Ref. [55], the effects of WDM suppression 
of the linear power spectrum and a suppression of the 
small-mass slope of the halo mass function were incorpo- 
rated into an estimate of the sensitivity of future weak 
lensing surveys to WDM. In Ref. [53], the WDM effects 
of the lack of small-mass halos in the mass function, po- 
tential cores in the halo density profile, and biasing of 
the smooth component were also included in estimates 
of the sensitivity to WDM in weak lensing surveys. 

In this paper, we incorporate results from our numer- 
ical simulations of WDM structure formation as well as 
the leading results from CDM halo statistics to construct 
a complete halo model of large scale structure in the case 
of WDM. We include new numerical simulations' mea- 
sures of the halo mass function, halo profile, and sub-halo 
mass function in this work. Importantly, since WDM is 
known to strongly affect halo substructure, we include 
the effects of halo substructure suppression here, which 
has not been included in prior work. We also incorpo- 
rate the change of the halo-profile density slope and the 
halo-profile concentration relation, which are known to 
depend on the initial matter power spectrum 2 1 [55] , as 
well an estimate of the biasing of the smooth dark mat- 
ter component in WDM similar to that in Ref. [53]. No- 
tably, we do not include cored halos in the case of WDM 
since these have not been observed in our simulations or 
others, nor are significant cores expected from analytic 
Gaussian-peak statistics methods [5j3 [57] . 

Two primary candidates for WDM particles are grav- 
itinos [55] and sterile neutrinos [55]. For concreteness, 
we employ the popular Dodelson-Widrow (DW) [29] 
thermally-produced sterile neutrino dark matter as our 
WDM particle dark matter candidate. All particle dark 
matter masses to which we refer are that of the DW ster- 
ile neutrino, unless stated otherwise. 



II. SIMULATIONS 



initial conditions include particle velocities due to the 
gravitational potential using the Zeldovich approxima- 
tion but we do not add random thermal velocities ap- 
propriate for WDM to the simulation particles. Gravi- 
tationally bound dark matter halos arc identified using 
the AMIGA'S Halo Finder software [33]. Each simula- 
tion set consists of a single realization of the density field 
but with varied power spectra of fluctuations appropriate 
for cold dark matter (CDM) and WDM cosmologies. We 
use the WDM transfer function in Ref. [7J valid for grav- 
itino particles. In Ref. [33] the transfer function of DW 
sterile neutrinos is shown to be nearly identical to that 
for gravitinos but with a scaling relationship between the 
particles masses. We characterize our WDM simulations 
in terms of the DW sterile neutrino mass. 

To investigate the halo mass function our first simula- 
tion set consists of 512 3 particles in a 60 3 Mpc 3 comoving 
box with mass resolution 5.7 x 1O 7 M and force resolu- 
tion 1.2 kpc. We run simulations for CDM and WDM 
cosmologies with sterile neutrino particle masses of 0.2, 
0.7, 1.7, and 4.4 keV. 

We use the high resolution simulations described in 
Ref. [T5] to investigate the subhalo mass function. These 
simulations use a 'zoom' technique to sample small vol- 
umes with fine mass resolution containing two Milky 
Way-sized halos, M ~ 2 x 10 12 M o , embedded in coarsely 
sampled larger volumes. The coarse volume is a 90 3 Mpc 3 
comoving box with ~ 510 3 particles with mass and force 
resolutions 9.2 x 10 4 Af Q and 275 pc in the refinement 
region. These simulations include CDM and 4.4, 11, and 
28 keV sterile-neutrino particle mass WDM cosmologies. 

The purpose of our final simulation set is to investigate 
the halo density profile. No differences were observed 
in the density profile slopes of the Milky Way-sized ha- 
los simulations, however these halos are well above the 
scale of power suppression in the WDM cosmologies. To 
reduce the halo mass to ~ 10 8 -M Q we selected one of 
these halos for resimulation in a smaller box, 4.5 3 Mpc 3 , 
with ~ 255 3 particles with mass and force resolutions of 
92Mq and 14 pc in the refinement region. We perform 
simulations for CDM and 28, 48, and 70 keV standard 
sterile-neutrino particle WDM cosmologies. We examine 
the halo density profiles at redshift 1.08 where the halo 
appears to be well relaxed in all cosmologies. 



We employ three sets of simulations using cosmolog- 
ical parameters consistent with the third year WMAP 
release [3D], spectral index n — 0.951, matter density 
n m — 0.238, baryon density f2& = 0.04, neutrino density 
f2„ = 0, Hubble parameter h = 0.73, and mass fluctu- 
ation with R = 8/i -1 Mpc, erg = 0.751. We assume all 
matter in the simulations is dark matter only but we use 
Ob for calculating the effects of baryons on the matter 
power spectrum. The simulations are conducted with 
the TV-body code GADGET-2 [3T] with initial conditions 
generated with the GRAFIC2 software package [35]. Our 



III. MATTER POWER SPECTRUM 

A. Cold Dark Matter 

The cold dark matter (CDM) power spectrum has been 
well-defined in the halo model of large scale structure 
(e.g., Ref. [15)). The dimensionless power spectrum is 

A 2 (fc)^^P, (3.1) 
2tt z 
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and the mass variance within region R is a 2 (R), 
dk 



k 



A (k) W 2 (kR) , 



(3.2) 



where W (kR) is the Fourier transform of the top-hat 
function with radius R. The quantity a is also given in 
terms of the halo mass given the mass enclosed in radius 
R for the mean density of the universe. 

We take the linear matter power spectrum for cold 
dark matter from Ref. [35]. The power spectrum is the 
primordial power spectrum modified by the appropriate 
transfer function, 



A 2 {k,z) = 5l 



ck 
Ho 



3+n 



T 2 (k,z) 



(3.3) 



where T(k,z) is the CDM transfer function, 8h is the 
amplitude of perturbations on the scale of today's horizon 
(and is absorbed into the normalization constant), and n 
is the initial power spectrum index. We employ the trans- 
fer function for CDM plus baryons from Ref. 35 1. This 
transfer function is sufficiently accurate for our purposes. 
The cosmological parameters we use for the transfer func- 
tions and analytic halo model are: spectral index n = 1, 
CMB temperature X"cmb = 2.726, dark matter density 
^dm = 0.23, baryon density — 0.04, neutrino density 
f2„ = 0, Hubble parameter h = 0.7, and mass fluctua- 
tion with R — 8 h^ 1 Mpc, as — 0.8, which are consistent 
with WMAP7 PP. We explore results at redshift z = 0. 



B. Warm Dark Matter 

For the WDM case, for concreteness, we use the ster- 
ile neutrino transfer function from Ref. [36 . There is 
a simple scaling relation to connect between sterile neu- 
trinos and gravitino WDM transfer functions given in 
Ref. [36] . The sterile neutrinos have an initial velocity 
dispersion, which allows them to escape the gravitational 
well of the small scale perturbations upon entering the 
horizon. This suppresses the linear matter power spec- 
trum for small scales (large A;), as shown in Fig. [I] The 
scale and magnitude of the suppression increases with 
velocity and, therefore, decreases with the mass of the 
sterile neutrino. The linear matter power spectrum for 
the suppressed sterile neutrino case at small scales is re- 
lated to that in CDM by: 



T s (fc) 



' -Psterile (&) 
-PCDM (fc) ' 



(3.4) 



In Ref. [36j , the linear matter power spectrum for ster- 
ile neutrinos was calculated using cambQ using a non- 
thermal momentum distribution and the evolution equa- 
tions for massive neutrinos. The resulting fit to the sterile 
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FIG. 1. Shown is the linear matter power spectrum as a func- 
tion of wavenumber for the CDM model and several WDM 
particle masses, as listed in the legend. The smaller the WDM 
particle mass, the greater the suppression of the power spec- 
trum on small scales (large k). 



neutrino transfer function is 

T s (k) = (1 + (afc)T M ■ 



(3.5) 



where v = 2.52, p = 3.08, and a is a function of the 
sterile neutrino mass, 
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a = a \ — 
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where a = 0.188, b 



h- 1 Mpc, (3.6) 



26/ \0.7 y 
0.858, c = -0.136, and d = 0.692. 



C. Nonlinear Matter Power Spectrum 

The nonlinear matter power spectrum consists of two 
parts: the one halo term and the two halo term [15j . 
These are denoted as and P 2 h, respectively: 



P (k) = P lh (fc) + P 2h (k) 



(3.7) 



The one halo term is Fourier transform of the two-point 
correlation function for two points that are inside the 
same halo: 



P lh (fc) 



dM- 



dn 
dM 



\u(k | M)| z 



(3.8) 



The one halo term depends only on the halo mass func- 
tion and the halo density profile. The two halo term is 
the Fourier transform of the two-point correlation func- 
tion for points that are in different halos: 

dn Mi 



P 2h (k) 



dMx 



dM- 



dM x p 
dn M 2 



■dM 2 p 
xPh h (k | M U M 2 



u(k | Mi) 
u (k | M 2 ) 



(3.9) 
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FIG. 2. Shown is the difference between the WDM models 
and CDM for the nonlinear matter power spectrum with only 
a change in the linear matter transfer function. The WDM 
particle masses are listed in the legend. This results in sup- 
pression of the power spectrum on small scales (large k) in 
the WDM models. 



In addition to dependence on the mass and density func- 
tions, the two halo term also depends on the power spec- 
trum of halos with masses Mi and Mi (Phh) ■ Phh can 
be approximated by the bias for each halo and the linear 
matter power spectrum: 

P hh (fc | M U M 2 ) « h (Mi) b 2 (Af 2 ) Pu n (k) . (3.10) 

In the nonlinear matter power spectrum, the suppres- 
sion of the linear matter power spectrum only directly 
effects the 2-halo term. As with the linear matter power 
spectrum, the nonlinear spectrum is suppressed at small 
scales, as shown in Fig. [2j 



IV. HALO BIAS FUNCTION 

The halo bias function describes how halos cluster rel- 
ative to the matter power spectrum. It is defined as 
the ratio of the halo power spectrum to the linear power 
spectrum [37], 



b 2 (*) 



Phh (fc) 

Pin W 



(4.1) 



As shown in Fig. [3j small mass halos are less strongly 
biased than higher mass halos. Ref. [35] showed that 
halos with M < M+ are more strongly clustered and halos 
with M > M+ are less strongly clustered than given from 
the bias function based on the standard Press-Schechter 
formalism 39 . M* is the typical mass scale of halos that 
are currently collapsing. Ref. [ID] also studied the halo 
bias function. They find that previous forms of the bias 
function overestimate the bias for halo masses slightly 
less than M+ and that the bias is approximately constant 
for small mass halos (M < 0.1M*). 
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FIG. 3. Shown is the halo bias as a function of mass for CDM 
and several WDM particle masses, as given in the legend. The 
difference in the WDM models comes only from the change 
in a (M). 



The halo bias function using a larger suite of simula- 
tions consistent with extended Press-Schechter formalism 
was found in Ref. [37], which we employ here. Specifi- 
cally, we use the bias function in §3.1, Eq. (6): 



b{v) = 1 — A- 



Bv h + Cv c 



(4.2) 



with parameters given in Ref. [37], Table 2, where v is 
defined by 

(4.3) 



cr(M) 



where a (M) is the variance at a given mass scale. 

The quantities r5 c and M* are defined in the spherical 
collapse model (see, e.g. Ref. [H]). Here, M* (z) is de- 
fined as the typical mass that is collapsing at redshift z. 
The value for M* is given by 



a(M it (z)) = 



D( Z y 



(4.4) 



where D is the linear growth, and note that D (0) = 1. 
The constant 6 C is the density for collapse in the spherical 
collapse model, and is S c — 1.69. 

In the case of WDM, M+ is not well-defined for all 
particle masses. For light masses, a (M) is less than <5 C 
for all masses (see Fig. [4]). This means that there are 
no typical halos that are currently collapsing. Since that 
is certainly not the case in the universe, these particle 
masses of WDM are not realistic. 

The halo bias calculated from a (M) for the CDM 
model and several WDM particle masses is shown in 
Fig. [3] The bias is defined in terms of v which, in 
turn, depends on a(M). Therefore, changes to the lin- 
ear matter power spectrum also affect the bias function. 
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find similar but more accurate relations for the halo mass 
function. We use the halo mass function from Ref. [41) . 
Appendix C: 



dn / \ Po dhia 
dM ^ 9 ^> M d~M 

with g (a) defined as 

,(.) = 5 ((f)" d + ^)exp(-| 



(5.3) 



(5.4) 



Their results for the constants B, d, e, / and g at sev- 
eral values of the mean interior density A are given in 
Ref. [H], Appendix C. 



FIG. 4. Shown is a (M) for the CDM model and several 
WDM particle masses, as given in the legend. Here, M sat- 
isfying a (Af ) = 5 C is the typical mass of halos which are 
currently collapsing. As shown, for very small WDM particle 
masses, such as 0.1 keV, typical fluctuation-scale halos have 
not collapsed, and little structure would have formed today. 



Since the halo bias function has not been studied in 
WDM structure simulations, we use the bias function 
from Ref. j37] for both the CDM and WDM models. 
Though the halo bias should be tested in the case of 
WDM with simulations, it is beyond the scope of the 
present work. 



B. Warm Dark Matter 

We find the halo mass function for CDM and three 
WDM particle masses in our simulations, described in 
pl| Our data shows that the suppression of the power 
spectrum on small scales, as expected, results in fewer low 
mass halos. We removed the upturn in number of halos 
with the smallest masses, which is a numerical artifact 
due to the discrete nature of the simulation We ht 
the WDM mass function as 



dn' w 
dM 



= 1 



M f 
M 



~ v dnc 
dM' 



(5.5) 



V. HALO MASS FUNCTION 

A. Cold Dark Matter 

The halo mass function, dn/dM is the comoving num- 
ber density of halos as a function of halo mass M and 
redshift z. From the spherical collapse model, overden- 
sities in the matter density field must be at or above 
the critical density S c to collapse. (Recall the definition 
of M* in equation (4.4 1.) The halo mass function can 
be estimated from the number of overdensities above the 
critical density. The following form was developed for the 
mass function, given Gaussian fluctuations in the initial 
density field: 



m 2 dn dM 

j;dM {M ' z) w 



dv 



(5.1) 



The variable v is defined in Eq. (4.3 1, and p z is the av- 



erage density of matter at redshift z. The function / {v) 
is defined by 



vf (v) 



exp 



(5.2) 



The PS halo mass function is an approximation. Moti- 
vated by extensions to PS theory, numerical simulations 



which does not include the erasure of the smallest scale 
halos. Here, dnc/dM is the CDM mass function, and Mf 
is the filtering mass, which is defined as 



M f 



47T _ 7T 3 



(5.6) 



where kt is the wavenumber at which the WDM transfer 



function, Eq. (3.5), has an amplitude of 1/2. For sterile 
neutrinos with masses 0.5 keV, 1 keV, 3 keV and 10 keV, 
kf ss 1.6, 2.8, 7.3, and 20.5 /iMpc -1 , respectively, and the 
filtering masses are 2.1 x 10 12 h~ 1 M Q , 3.6 x 10 11 h~ 1 M Q , 
2.1 x 10 10 /i _1 M Q and 9.8 x 10 8 h~ 1 M G respectively. 



Eq. (5.5) is an accurate fit to the simulation's mass func- 
tion for 77 = 1.2 (see Fig. [5|. 

Since the WDM particles free-stream out of the 
small scale perturbations, dark matter halo formation is 
strongly suppressed below the free-streaming scale i?f s , 
which is defined as 



R [s = 0.2 (fi dm /i 5 



,1/3 _4/ 3 



(5.7) 



This corresponds to a free-streaming scale mass Mf s at 
mean density (not to be confused with the filtering mass 
Mf). We use this mass for our small scale cut-off in 
the mass function. We use an expression for the small 
halo suppression effect that is approximately unity when 
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FIG. 5. Shown is the halo mass function from our simulations 
and our fits to that data. The thick lines are the data from 
each simulation, as listed in the legend. The thin lines are the 
CDM data multiplied by our fitting factor, which is given in 
Eq. (pT5||. 
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FIG. 6. Shown is the halo mass function for the CDM model 
and several WDM particle masses, as given in the legend. The 
mass functions are normalized to the CDM mass functions at 
large masses. Our WDM mass functions include a suppression 
for small M and a small scale cut-off. 



in 3> Mf s and goes to zero when m <C Mf s to provide a 
continuous cut-off. Our WDM mass function is then 



Ilm 



M 2 



M 2 



Ml 



~M 



dn c 
dM' 



(5.8) 



The full mass function for our CDM and WDM models 
is shown in Fig. [6] The effect of just the change to the 
mass function on the nonlinear matter power spectrum is 
shown in Fig. [7} The decrease in small mass halos results 
in the suppression of the power spectrum at small scales. 



FIG. 7. The percent difference between the nonlinear matter 
power spectrum for the CDM model and several WDM parti- 
cle masses, as listed in the legend, is plotted versus wavenum- 
ber k. In this graph, the only difference between the CDM 
and WD M m odels is the change to the mass function, as given 
in Eq.s (5.81. The decrease in small mass halos results in the 



suppression of the power spectrum at small scales. 

VI. HALO DENSITY PROFILE 
A. Cold Dark Matter 

A general halo density profile is given by Ref. 
p{r\M) = 



G8-7)/V 



(6.1) 



For CDM, we employ the Navarro, Frenk and White 
(NFW) profile 01], which has a = 1, /3 = 3 and 7 = 1. 
The radius r s is the radius at which the density function 
has a logarithmic slope of —2. For a spherically symmet- 
ric density profile, the mass of the halo is 



M = I drAirr 2 p (r | M) , 
/o 



(6.2) 



for the virial radius i? v ir- For the NFW profile, the in- 
tegral has an analytic solution, which specifies the value 
of p s : 



M = 47rp s rf In (1 + c) 



c 

l + c 



(6.3) 



The concentration c is defined as c = R v - lr /r s . When 
calculating the nonlinear matter power spectrum, we use 
the Fourier transform of the density profile, 



u(k\M) 



Rv " , , 2 sin(fcr) p(r I M) 

dr4nr 2 — ^ — '- ' y J. ' . (6.4) 
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B. Warm Dark Matter 

The density profile of halos in WDM simulations has 
previously been studied by Avila-Reese and Colin ct 
al. O US H6]. In [45], Milky Way size halos are sim- 
ulated for WDM particles of mass 0.6, 1. and 2. keV. 
Ref. [35] found that the density profile of the halos is de- 
scribed by the NFW profile, but WDM halos can have an 
inner slope slightly shallower than — 1. In Ref. 46], they 
simulate CDM and WDM halos with masses down to 
O.OlAff. These halos were found to also be well described 
by the NFW profile. For WDM halos with masses below 
Mf, it was found that the inner slope is shallower than in 
comparable CDM halos. In Ref. [33], halos with masses 
close to the filtering mass Mf were simulated, which found 
that the innermost (r < 0.02i? v ; r ) logarithmic slope of 
WDM halos is steeper than in the CDM model. Outside 
of this volume, the density profile for WDM halos is shal- 
lower than the NFW fit. Ref. gSJ finds that, using the 
density profile of Eq. (6.1 1, a cored profile with a = 0.7, 
j3 = 3 and 7 = fits the WDM halos. Ref. [53] uses a 
halo density profile where the size of the core increases 
as the mass of the halo decreases. The small mass halos 
(M < 10 11 h^ 1 M ), in Ref. [23] have cores with radii on 
the order of r s . On the other hand, more recent work in 
Refs. [25J [2T] studies whether WDM halos can have such 
cores, and shows that, in halo collapse modeling, the core 
of a WDM halo is smaller than r < 10~ 3 i? vir . Below this 
scale, a core may need to be included in the density pro- 
file. Therefore, it is unlikely that WDM halos have large 
cores. 

WDM halos that are smaller than the filtering mass 
(see Eq. ( |5.6| )) have a shallower inner slope compared 
to CDM halos gSJ gB]. Ref. [25] shows that, in CDM 
simulations, the slope of the inner profile depends on the 
effective spectral index n e ff of the initial power spectrum 
of the density perturbations at the scale of the power 
spectrum sampled by the halo fchaio- That is, P (fchaio) oc 

fc hafo- Ref - E3 flnds that 



7 = 



9 + 3n eff 



(6.5) 



where 7 is an exponent from Eq. (6.1 ). This is relevant to 



our WDM model since the suppression of the small-scale 
structure acts to decrease the effective spectral index at 
small scales, as in Fig. [I] This is what leads to a decrease 
in 7, the inner slope of the profile, though since n e s is not 
well defined in WDM for a fixed halo mass, one cannot 
specify 7 analytically. 

We use our simulation of a 10 8 Mq halo for CDM and 
three sterile neutrino WDM particle masses: 28 keV, 48 
keV, and 70 keV. The density profile of these halos are 
shown in Fig. [8] We find that the simulated CDM and 70 
keV WDM halos have an inner slope of —1.2. The inner 
slope of the simulated 48 keV WDM halo is —1.0, and the 
simulated 28 keV WDM halo has an inner slope of —0.8. 
We use an inner slope of —1 for the CDM cosmology, 
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FIG. 8. This is a plot of the density versus radius for a 10 s Mq 
halo simulated in four different WDM cosmologies, as given in 
the legend. The thick lines are data from the simulation, and 
thin lines show the estimated inner slope for the halo density 
profile, with 7 = 1.2, 7 = 1.2, 7 = 1.0, 7 = 0.8 for the CDM, 
70 keV, 48 keV and 28 keV cases, respectively. The value 7 
describes the profile inner slope in Eq. (16. 11. 



as in the NFW profile. To model the effect of the sup- 
pression of the density profile, the inner slopes of halos 
composed of 70 keV, 48 keV, and 28 keV WDM particles 
are chosen to be —1, —0.8 and —0.6, respectively. We 
use these values to interpolate the inner slope for other 
WDM particle masses. 

To include this effect, we make changes to the expo- 
nents in the profile density function (Eq. (6.1)). The 
exponent 7 describes the logarithmic slope of the cen- 
tral section of the halo, where r/r s <C 1. For the inner 
slope to be shallower, we use 7 < 1. The smaller the 
mass of the WDM particle, the shallower the inner slope. 
When the inner slope of the density profile decreases, the 
density of the central region of the halo increases for a 
constant halo mass. 

The exponent (3 describes the slope of the density pro- 
file for the outer edge of the halo, where r/r s 1. The 
quantity a describes the sharpness of the change between 
the inner and outer sections of the density profile. For 
larger a, the changeover between a slope of 7 and a slope 
of (3 is sharper. It takes place within a smaller region 
around r/r s = 1. 

To maintain the definition of r s as the radius at which 
the density profile has a logarithmic slope of —2, the 
exponents are constrained to be: 



1 = 



7 



2-/3 



= 4-7. 



(6.6) 



Note that this relation is valid for a = 1, j3 = 3 and 
7=1, which is the NFW profile. 

Since the profile exponent change in 7 changes the den- 
sity distribution, in order to maintain a constant halo 
mass, the value of p s must be rescaled with the pro- 



8 



~i 1 — n~| 1 1 — n~| 1 1 — rr~| 1 1 — rr~| 1 1 — r 



10° 



10-' 



io- : 



io- 3 



CDM 
10 keV 
3 keV 
1 keV 
0.5 keV 



10' 



10 s 



10 9 



10 10 10 1 
m [h 1 M D ] 



10 1; 




10° 10 1 
k [h Mpc 1 ] 



FIG. 9. The Fourier transform of the halo density profile 
u (k\M) for a constant wavenumber of 10 2 /i _1 Mpc is plotted 
as a function of halo mass M for the models listed in the 
legend. 



file changes, which we incorporate in the WDM models. 
Shallower 7 requires an increase in p s , which broadens 
the density profile for a fixed halo mass. The Fourier 
transform of the density profile is employed in the non- 
linear power spectrum. Fig. [9] is a plot of u (k\M) for a 
constant wavenumber k. The increase in density around 
the halo center increases u (k\M) for large k given a con- 
stant M. Fig. 10 shows the effect of changes to the halo 
density profile on the nonlinear matter power spectrum. 
The broadening of the central part of the halo and sup- 
pression of the central peak profile results in a transfer 
of power from the smallest scales to slightly larger scales, 
seen in Fig. [l0| 



VII. CONCENTRATION 



The concentration of an NFW (see \ VI I CDM halo has 
been found to be [24] 



c (to, z) 



1 + Z \771* (z) 



-0.13 



(7.1) 



We apply the methodology described in Ref. [23] to the 
WDM linear power spectrum to determine the change in 
halo concentration. 

The model in Ref. |24] is determined by two equations 
and the parameters F and K: 



M* (z c ) = FM, 



and 



c (M, z c 



K 



l + z c 



1 



(7.2) 



(7.3) 



FIG. 10. The percent difference between the nonlinear matter 
power spectrum for the CDM model and several WDM parti- 
cle masses, as listed in the legend, is plotted versus wavenum- 
ber k. In this graph, the only difference between the CDM 
and WDM models is the change to the inner slope of the halo 
density profile. The broadening of the central part of the halo 
and suppression of the peak results in a transfer of power from 
the smallest scales to slightly larger scales, as shown by the 
peak in the above plot. 



We assume that a (M) is approximately a power law 
at a mass of FM+: 



c(xM~ a ". 

Then, we can calculate a c from a (M) : 

_ d (in (j ^) _ —M da 
" c ~ d (In M) ~ ~d~M' 



(7.4) 



(7.5) 



For the CDM case, the best fit parameters for F andif 
are 0.01 and 4.0, respectively. Then, a c — 0.13, as in 
Eq. j7\I]). This is valid for o'.OlM* < M < 100M* [H]. 
11 shows a r 



Fig. 



as a function of WDM particle mass. 



where M± is defined in Eq. (4.4) 



The resulting concentration for halos in the CDM models 
and for several WDM masses is shown in Fig. [12J Since 
the exponent a c decreases with WDM particle mass (see 
Fig. [TTj ), the dependence of the concentration on halo 
mass also decreases. In addition, the concentration of 
halos in WDM models is less than in CDM [H SB US- 
Since a c depends on the linear matter power spec- 
trum through cr(M), changing the concentration has 
an effect only if the linear matter power spectrum also 
changes. Fig. [F3| shows the effect of the WDM linear mat- 
ter power spectrum and concentration on the nonlinear 
matter power spectrum. The decrease in the concen- 
tration (see Fig. 12) increases the scale radius r s , from 
Eq. (6.1 ), which stretches the inner part of the halo pro- 
file relative to the radius. This also results in a transfer 
of power from the smallest scales to slightly larger scales, 
as shown by the peak in the nonlinear matter power spec- 
trum. 
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FIG. 11. The exponent in the concentration equation, 
Eq. (7.4|, is plotted as a function of the WDM particle mass. 
As particle mass increases, a c approaches 0.13, its value in 
the CDM model. 
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FIG. 12. The concentration of halos is plotted as a function 
of halo mass for the CDM model and several WDM parti- 
cle masses, as given in the legend. Since the exponent a c 
decreases with WDM particle mass (see Fig. Ill, the depen- 
dence of the concentration on halo mass also decreases. 



VIII. SUBSTRUCTURE 

Clearly, in WDM, small mass halos are suppressed in 
number, and therefore subhalos will be suppressed as 
they are accreted halos in the hierarchical formation of 
large scale structure. Dolney et al. |17j study the effects 
of substructure on the CDM halo model. A halo con- 
sists of a smooth mass component and several subhalos. 
Most of the mass of the halo is in the smooth component. 
A fraction / of the mass is in the subhalos. Ref. [47] 
uses / ~ 0.1 and assumes that all subhalos are smaller 
than 0.01M where M is the mass of the host halo. With 
substructure, the one halo term becomes a sum of four 
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FIG. 13. The percent difference between the nonlinear matter 
power spectrum for the CDM model and several WDM parti- 
cle masses, as listed in the legend, is plotted versus wavenum- 
ber k. In this graph, the WDM models differ from the CDM 
model in both the transfer function and concentration. The 



decrease in the concentration (see Fig. 12 1 increases the scale 
radius r s , from Eq. (6.1 1, which stretches the inner part of the 
halo profile relative to the radius. This results in a transfer 
of power from the smallest scales to slightly larger scales, as 
shown by the peak in the above plot. 



different terms: 



-Pih — Pss + Psc + Pic + P2c 



(8.1) 



The first term in Eq. (8.1 1, P ss describes correlations be- 
tween two points of the smooth component of the host 
halo: 



P ss (k) 



dM 



.dN 
dM 



Ms 
? 



\U(k | M a )\* 



(8.2) 



Throughout this section, variables that are capital letters 
describe the host halo, and lower case variables are asso- 
ciated with the subhalos. Note the similarity of Eq 
to the one halo term without substructure, Eq. (13.8 1 



1) 



The quantity M s is the mass of the host halo that is in 
the smooth component {M s = M (1 — /)), dN /dM is the 
mass function for host halos, and U is the den sity profile 
for host halos. The second term in Eq. ( |8.1| ) describes 
correlations between the smooth component of the host 
halo and the subhalos: 

/* 7 TV T 71 r 

P sc (k) = 2 dM—^U (k | M s ) U c (k | M s ) 
J dM p 

x / dm^-u(k | m). (8.3) 
J dm p 

The expression U c is the density profile for subhalos 
within the host halo, and dn/dm and u are the mass func- 
tion and density profile for the subhalos. The third term 
from Eq. (8.1 1 is the correlation for two points within the 



10 



same subhalo: 



Pic (k) = J dM 



dN 
dM 



dm p(™) (8.4) 
dm \ p 



The fourth term is the correlation between two different 
subhalos in the same host halo: 

P 2c (k) = J dM^\U c (k\M s )\ 2 



dm- —u (k I m) 

dm p 



(8.5) 



Dolney et al. [47] assume that the effect of substructure 
is negligible on the scale of correlations between two dif- 
ferent halos, and we do the same. Therefore, the two 



halo term (Eq. (3.9)) remains the same. 



The host halo mass function dN/dM and density pro- 
file U are as defined in the previous sections. For sub- 
structure, we need expressions for the spatial profile of 
subhalos within the host halo U c , the mass function for 
subhalos dn/dm and the density profile with the subhalos 
u. We use the same density profile for the subhalos as for 
the host halos. The distribution of subhalos in the host 
halo is taken to be the density profile of the host halo U. 
The concentration for the subhalos is given by [23] : 



c (m, z) = 



7.5 



m 



-0.30 



(8.6) 



The exponent —0.30 corresponds to F w 200 where F is 



taken from Eq. (7.2). The subhalo mass function is 



dn ,. , / M\' i dm 

- — (M) dm = No I — I — . 
am V m ) m 



(8.7) 



The constant Nq is determined by the fraction / of the 
halo mass in subhalos. This relation can be solved ana- 
lytically to find 



/ = 



, m dn O.Ol^-^iVo 

dm— — (M) = -. 

Mdm K ' 1-fj, 



(8.8) 



where the most massive subhalos are assumed to be one 
hundredth of the mass of the host halos. Dolney et al. [37J 
usea/iK 0.9 for CDM. 

A. Warm Dark Matter 

In the CDM model, we employ the same density pro- 
file for halos and subhalos. When considering the effect 
of WDM on the subhalo density profile, we alter it in 
the same was as the halo density profile, which is de- 
scribed in §VT| Fig. [14] shows how changing the subhalo 
density profile with WDM particle mass affects the non- 
linear matter power spectrum. In this figure, the CDM 
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FIG. 14. The percent difference between the nonlinear matter 
power spectrum for the CDM model with substructure and 
several WDM particle masses, as listed in the legend, is plot- 
ted versus wavenumber k. In this graph, the WDM models 
differ from the CDM model in only the subhalo density pro- 
file. As with changes to the halo density profile, decreasing 
the inner slope of the subhalo density profile decreases the 
power on the smallest scales (large k). 



model also includes substructure. Since the center of the 
subhalo clumps shallower cusps, the smallest scales (high- 
est k) transfer power to slightly larger scales. But, those 
scales also lose power to this effect, so the net result is 
negative. 

For the subhalo concentration, we follow the proce- 
dure outlined in §VII[ but use F s» 200 where F is taken 
from Eq. (7.2). This allows us to calculate the exponent 
a c , in Eq. (7.4), for the subhalo concentration. For the 



WDM particle masses 0.5 keV, 1 keV, 3 keV, and 10 keV, 
this exponent is 0.23, 0.25, 0.26, and 0.26, respectively. 
The effect of the subhalo concentration on the nonlinear 
matter power spectrum is shown in Fig. 15 The con- 



centration of the subhalos decreases, which decreases the 
amplitude of the power spectrum on the smallest scales 
(large k). 

For the subhalo mass function, we use our simulations 
to determine the change in the subhalo mass function 
relative to WDM particle mass. We want a function such 
that 



am 



dnc 
dm 



(8.9) 



where m! — m/Mf and Mf is the filtering mass. To find 
g(m'), we take the ratio of the WDM to CDM mass 
functions and change variables to m! . The result of this 
procedure for both the halo and subhalo mass functions 
is shown in Fig. [16] As expected from |VB[ the halo 
mass function matches the multiplicative factor given in 
Eq. (pTBl That is, 



dN w /dM 
dNc/dM 



(1 + M'- 1 )- 



(8.10) 
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FIG. 15. The percent difference between the nonlinear matter 
power spectrum for the CDM model with substructure and 
several WDM particle masses, as listed in the legend, is plot- 
ted versus wavenumber k. In this graph, the WDM models 
differ from the CDM model in the transfer function and the 
subhalo concentration. The concentration of the subhalos de- 
creases, which decreases the amplitude of the power spectrum 
on the smallest scales (large k). 



where M' = M/Mf of the host halo. When the same 
procedure is applied to the subhalo mass function, we find 
that the data fits the same expression as for the halo mass 
function. As with the host halos, we use a small scale cut- 
off for the subhalos. Our subhalo mass function changes 
with WDM particle mass i n exa ctly the same wa y as the 
halo mass function in Eq. (5.8), but using Eq. (8.7) for 



the CDM subhalo mass function. There may be slight 
evidence, as seen in Fig. |16[ that the mass function for 
subhalos is less steep as a function of m! . 

We normalize the subhalo mass function so that there 
are the same number of the largest subhalos (0.01 the 
halo mass) in CDM and WDM. The constant An from 



Eq. ( |8.7[ ) in the WDM model is related to CDM model 
by 



AW = I 1 + 100 



.M^ r < 
' M 



0c ■ 



(8.11) 



The fraction of the halo mass in subhalos is determined 

) . Note the dependence on the 
ler the mass of the host halo, 



by the integral in Eq. (8.8 
host halo mass 



The sma 

the smaller the fraction of the mass in subhalos. 

The effect of the subhalo mass function on the nonlin- 



ear matter power spectrum is shown in Fig. 17 As the 
WDM particle mass decreases, the slope of the subhalo 
mass function p decreases, and Mf s increases. This de- 
creases the fraction of the host halo mass in subhalos and 
increases the fraction of the mass in the smooth compo- 
nent of the host halo. This is the reason for the increase 
in power in the nonlinear matter power spectrum between 
0.1 and 1 h Mpc -1 . Power is transferred from the small- 



est scales, the subhalos, to somewhat larger scales, the 
host halo. 



IX. SMOOTH BACKGROUND 

The suppression of small mass halos results in some of 
the dark matter being in a smooth background compo- 
nent instead of collapsing into halos. The fraction of the 
dark matter in halos is: 



1 

Po 



dMM 



dn 
dM' 



(9.1) 



The fraction of the mass in halos is 0.95 for a 10 keV 
WDM particle, 0.76 for a 3 keV particle, 0.52 for a 1 keV 
particle and 0.36 for a 0.5 keV particle. We assume that 
the smooth component of the dark matter is related to 
the linear matter power spectrum by 



,2 _ -Pss (k) 

where b s is derived in Ref. [23 : 

1 - fhb eS 



6 S = 



1 - A 



and b eS is the effective mass weighted halo bias: 

e _ J dMM (dn/dM) b (M) 
/ dMM (dn/dM) 



(9.2) 



(9.3) 



(9.4) 



The power spectrum for two points in the smooth com- 
ponent is [33J 



P ss (fc) = b 2 s P lin (k) . 



(9.5) 



and the correlation between the smooth component and 
halos is 



Psh (k) 



b s Pun (k) 



Ph 



uTI 

dM M u(k I M) b(M), 

dM y 1 ; v ' 



(9.6) 



where ph = fhP- Then, the full power spectrum is: 
P (k) = (1 - f h ) 2 P ss (k) 

+2 (1 - / h ) APh (k) 
+Pi h (k)+P 2h (k). (9.7) 

The full power spectrum and its components for a 0.7 



keV WDM particle are shown in Fig. 18 Note that the 
smooth-smooth and smooth-halo terms are much smaller 
than the 2-halo term. Therefore, these terms should not 
significantly affect our results. 



X. RESULTS 

Fig. [19] is a graph of the percent difference of the non- 
linear matter power spectrum in our WDM model with 
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FIG. 16. The ratio of the WDM to CDM halo and subhalo mass functions from our set of simulations is plotted as a function 
of the halo or subhalo mass m' = m/Mf for several WDM particle masses, as given in the legends. The ratio of mass functions 
as a function of m' fits the same form, listed in the legend, for both the halo and subhalo mass functions. 




FIG. 17. The percent difference between the nonlinear matter 
power spectrum for the CDM model with substructure and 
several WDM particle masses, as listed in the legend, is plot- 
ted versus wavenumber k. In this graph, the WDM models 
differ from the CDM model in only the subhalo mass func- 
tion. The changes to the subhalo mass function decrease the 
fraction of the host halo mass in subhalos and increase the 
fraction of the mass in the smooth component of the host 
halo. Power is transferred from the smallest scales, the sub- 
halos, to somewhat larger scales, the host halo. 

all effects except substructure. That is, the transfer func- 
tion, mass function, halo density profile, and concentra- 
tion are altered as described in the previous sections. At 
k s» lOh Mpc -1 , the strongest effects are from the con- 
centration and halo density profile, which increase the 
power at these scales. All four of our alterations to the 



FIG. 18. The power spectrum for a 0.7 keV particle mass, is 
plotted versus wavenumber k. This graph includes the WDM 
mass function and the effect of a smooth background com- 
ponent to the dark matter. The smooth-smooth, smooth- 
halo, 1-halo and 2-halo power spectrums for a 0.7 keV WDM 
particle are also shown. Note that the smooth-smooth and 
smooth-halo terms are much smaller than the 2-halo term. 



WDM model in this graph contribute to the decrease in 
power at small scales (large k). 

Fig. [20] is a graph of the percent difference of the non- 
linear matter power spectrum in our WDM model, in- 
cluding only the effects of WDM on the substructure and 
the transfer function, compared to the CDM model with 
substructure. The dominant effect at the 10 h Mpc -1 
scale, especially for the 10 keV WDM particle, is from 
the subhalo mass function. For the 0.5 keV WDM par- 
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FIG. 19. The percent difference between the nonlinear matter 
power spectrum for the CDM model and several WDM parti- 
cle masses, as listed in the legend, is plotted versus wavenum- 
ber k. This graph shows the results of changing the transfer 
function, mass function, halo density profile, and concentra- 
tion. It does not include substructure. At k « lOh Mpc - , 
the strongest effects are from the concentration and halo den- 
sity profile, which increase the power at these scales. All four 
of our alterations to the WDM model in this graph contribute 
to the decrease in power at small scales (large k). 



FIG. 21. The percent difference between the nonlinear mat- 
ter power spectrum for the CDM model with substructure 
and several WDM particle masses, as listed in the legend, is 
plotted versus wavenumber k. This graph includes all of the 
adjustments to the WDM model detailed in this paper. This 
figure greatly resembles Fig. |19| with the increase in power 
at k ~ lO/i Mpc -1 . Therefore, the largest effects are from 
the changes to the host halo properties, with changes to the 
substructure and the addition of a smooth background com- 
ponent to the dark matter being sub-dominant. 
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FIG. 20. The percent difference between the nonlinear mat- 
ter power spectrum for the CDM model with substructure 
and several WDM particle masses, as listed in the legend, is 
plotted versus wavenumber k. This graph shows the results 
of changing the transfer function, the subhalo mass function, 
subhalo density profile, and subhalo concentration. The dom- 
inant effect here, especially for the 10 keV WDM particle, is 
from the subhalo density profile. For the 0.5 keV WDM parti- 
cle, the largest effect is the decrease in power from the changes 
to the subhalo concentration. 



ticle, the largest effect is the decrease in power from the 
changes to the subhalo concentration. 

Fig. [21] is a graph of the percent difference of the non- 
linear matter power spectrum for our WDM model with 
all of the effects discussed in this paper compared to 
the CDM model with substructure. The features in the 
changes of the nonlinear matter power spectrum greatly 
resemble those in Fig. 
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Therefore, the largest effects 
are from the changes to the main halo in WDM, with 
changes to the substructure and the addition of a smooth 
background component to the dark matter being smaller 
effects. 



XI. CONCLUSION 

We have presented our work on a broad assessment 
of the effects of WDM in the halo model of large-scale 
structure in two-point statistics, employing results from 
our own set of simulations as well as previously published 
work. This framework may be incorporated in a number 
of applications, including forecasts for large scale struc- 
ture measures such as weak lensing [22, 23J, in galaxy 
clustering two-point function measures of the power spec- 
trum 48J. We have included the effects to the linear 
matter power spectrum, halo density profile, halo con- 
centration relation, halo mass function, subhalo density 
profile, subhalo mass function and biasing of the smooth 
dark matter component in the case of WDM. 

We find a drastic difference in the nonlinear matter 
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power spectrum predicted by the halo model between 
WDM and CDM even for reasonably "cold" WDM par- 
ticle masses (m s = 10 kcV) at the smallest scales, as 
expected. We show that, counter-intuitively, WDM pro- 
duces an enhancement of power at intermediate scales 
due to the softening of halo density concentrations and 
profiles. Since cored profiles are not found in WDM simu- 
lations, nor expected in analytic Gaussian peak statistics 
of WDM [2SJ[37], we do not include cored profiles, which 
differs our results from Ref. [23]. Also different from 
that work, we also include the effects of halo substruc- 
ture, subhalo density profiles and subhalo mass functions. 
Ref. [32] only included the effects of the halo profile and 
halo mass function, and assumed all mass to be in halos, 
which is not the case in WDM. 

We find that the host halo effects dominate the over- 
all effects of WDM versus CDM in the nonlinear regime. 
Though our detailed results differ from previous work, 
the overall magnitude of the effects in the intermediate 
scale (k ~ 10 — 100 h Mpc -1 ) regimes relevant for weak 
lensing are comparable, and therefore, the estimates for 



weak lensing sensitivity in that work are likely not dras- 
tically changed when incorporating our model in a weak 
lensing forecast. 

The halo model has broad applications to observable 
cosmological statistics. With the help of the nonlinear 
structure framework of WDM versus CDM like that pre- 
sented here, the observed large to small-scale clustering 
of matter, gas, and galaxies may be able to shed light on 
the nature of dark matter and its primordial origin. 
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